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1 . STATEMENT O F WORK 


The statement of work includes data analysis and supporting research in 
coimection with the following broad objectives: 

(1) Provide a precise and accurate geometric description of the 
earth's surface. 

(2) Provide a precise and accurate mathematical description of 
ihe earth's gravitational field. 

(3) Determine time variations of the geometry of the ocean sur- 
face, the solid earth, the gravity field and other geophysical 
parameters. 
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2. ACTIVITIES RELATED TO THE NGSP 
(Grant No. NGL 36-008-093) 

2. 1 Satellite Triangulation in Europe from WEST and ISAGEX Data 

In 1974 the Department of Geodetic Science at The Ohio State University 
(OSU) obtained observational data that was acquired during the Western European 
Satellite Triangulation (WEST) program and the International Satellite Geodesj'' 
Experiment (ISAGEX) campaign. 

The purpose of obtaining this data was twofold. Primarily, it was intended 
that a geometric solution be performed to improve the present values of coor- 
dinates of the European stations in the OSU WN14 solution. The secondary aim 
was to add some new stations and to assess the quality of the WN14 solution wth 
the help of the additional data available. 

2.11 Data 

The sets of data were thoroughly described in the Fifteenth and Sixteenth 
Semiannual Status Reports. No additional data has been acquired. 

2.12 Adjustments 

A total of three solutions were performed. Solution WEST 33 and WEST34 
contain only observational data of the Western European Satellite Triangulation 
and 1VEST-ISAGEX36 (W. 1. 36) is a combination solution containing also the data 
of the ISAGEX campaign, hi all cases only the single image data were processed. 
Preliminary computations with the seven image data (with assumed statistics 
about the observations) resulted in seemingly distorted station coordinates. In 
absence of any knowledge about the statistics of the observational data, it was 
felt very doubtful that the seven image data could improve the results of the 
single image data computation. Any further analysis with seven image data, 
tlierefore, must await loiowledge about the statistics of the data. 
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The results of the WEST33 solution were submitted the the XVI General 
Assembly of the lUGG held in Grenoble, France from August 18 through Sep- 
tember 6, 1975. The solutions WEST34 and W.I. 36 are documented in Reports 
of the Department of Geodetic Science No, 232 , OSU. 

The solutions WEST33 and WEST34 differ basically in the number of base 
lines which were constrained. It was recognized that the base line TROMSO- 
CATANIA was not sufficient to transfer scale to the whole net. The WEST 
satellite network is considered as consisting of two blocte: the central Euro- 
pean block with a large number of observations and the northern block which is 
connected to the central block by relatively few observations (namely, between 
TEOMSO and some stations of the central block) . An overall scale factor of 
10 ppm was computed between the ED 50 coordinates and the adjusted values . 
Comparing individual chords in the two systems, it became clear that all chords 
originating from TROMSO yield a significantly smaller scale factor. Also the 
scale for the central area is partly inherent in the weighted positional constraints 
of the WN14 stations. It thus became necessary to include more chord constraints, 
especially in the central area. These were taken from [Ehmsperger, 1974] . 

2. 13 Results 

(1) The ISAGEX data added three stations to the WEST34 system. Due to the 
small number of ISAGEX observations, only a minor improvement could be 
gained by the addition of the ISAGEX data. 

(2) The coordinates of about seven stations still exhibit extraordinarily 
large standard deviations. This is an immediate result of the increased 
variances of the observational data (see Table 4 in [Ehmsperger, 1974] for 
the standard errors stationwise, as obtained by the smootliing procedure), or 
it is due to lack of a sufficient nttmber of observations or to unfavorable geo- 
metrical conditions . 

(3a) Various transformations were carried out between the W. 1. 36 solution 
and other systems. Wliile the computed rotation angles between the W.I. 36 and 
the ED50 system are all below the one sigma level, the translation parameters 
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agree well with those of solutions of other investigators. 


(3b) Special effort was made to find the scale factor between the ED50 and the 
W.I. 36 system. The seven parameter transformation gives a scale factor of 
A(ppm) = 6.12± 2.77. However, there is evidence [Wei^itman, 19 75] that 
current publications give revised coordinates for the terminal station of the 
European base lines which are based on traverse adjustment. Deleting these 
stations and the stations with extraordinarily large variances in the coordinates 
(see 2), the scale factor was computed to A (ppm) = 5.92± 3.54. 

(4) The variances of the coordinates for the stations common to the ^VN14 
system were all significantly reduced. 


2.2 Similarity Transformation and Geodetic Nefvork Distortions Based on 

Doppler Observations 

The purpose of this investigation was twofold; (a) some theoretical con- 
tributions are given to the transformation models as used in geodesy, and (b) 
distortions in geodetic networks are investigated based on these transformation 
models . 

This investigation which is laid down in Reports of the Department of Geo - 
detic Science No. 235 , OSU (with the same title as this section), presents a review 
of die commonly used transformation models and provides a geometrical inter- 
pretation for most of these models . It is shown that the translation components 
as computed from the so-called Molodenskii Model for seven transformation 
parameters should not be interpreted as the vector between the origins of the two 
coordinate systems involved. Only the Bursa Model permits such an interpretation. 
It is furtlier shown that as for direct transformations, both the Bursa and Molo- 
denskii Models give identical results, i. e. , die same station coordinates and the 
same variance-covariance matrix. Also it is shown that the parameters as 
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computed from both models, differ only ia the translation components. An 
expression is given to' exactly convert one set of translation vectors into the other . 

Pitfalls in partial solutions are also discussed. It has been proposed in. 
several publications [Bursa, 1967 and Kumar, 1972] that the orientation angles 
and the scale factor be computed separately from direction and chord comparisons 
by using them in all combinations between stations. Based on the method of 
eliminating parameters, it is shown that these proposed approaches are wrong 
as far as the "selection" of direction and chords is concerned. Although these 
approaches yield parameters which are close to those computed without eliminating 
some of the seven parameters, their standard deviations are much too optimistic. 

It is shown that only as many directions and chords can be used in the computation 
as are needed to completely determine the shape or size of the polyhedron implied 
in the set of Cartesian coordinates . Each additional element (direction or chord) 
causes the normal matrix to be singular. 

In the closing section of this report a number of tables and maps indicating 
distortions in the North American Datum 1927, the Precise Traverse System 
M-R-72, the Australian Geodetic Datum and the South American Datum 1969 
are given. The residuals of the coordinates are scanned for systematic patterns 
after transforming the geodetic system to the NWL9D system. Also, an attempt 
has been made to plot maps showing scale distortions on the NAD27. 

Since the Reports of the Department of Geodetic Science No. 235 , OSU 
does not include tables of coordinates, aU the coordinates which were available 
for this investigation are given below. Some of these coordinates were already 
given in the Fifteenth Semiannual Status Report. The NWL-9D coordinates 
have been corrected in the meantime as explained in Attachment 1 of this report. 
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NAD 1927 Coordinates for North American Stations 
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20176 

AJCl, ARIZONA 

32 

26 

54.600 

247 

08 

51.637 

399.77 

20177 

DOUGLAS. ARIZONA 

31 

22 

36.952 

250 

27 

52.654 

1197.44 

20208 

KINGMAN. ARIZONA 

3 5 

11 

^-^8. 172 

245 

57 

31.854 

1 1 1 ? . 44 

30030 

GREEN RIVER, UTAH 

38 

58 

44.214 

249 

53 

17.276 

1288.256 

51004 

charlfston,hv 

38 

22 

30.835 

278 

24 

32.776 

249.38? 

51005 

CORBIN, KENTUCKY 

36 

57 

21.781 

275 

53 

08.520 

354.976 

51006 

CLEVELAND, TENNES. 

35 

09 

06.961 

275 

06 

55. 188 

263 .137 

51007 

LAURENS, south CAR. 

34 

35 

08.446 

277 

56 

35.479 

182,54.1 

51009 

SHELBY. ALABAMA 

33 

07 

03,355 

273 

30 

00.801 

181.752 

51010 

sandersville.g. 

33 

03 

38.462 

277 

05 

29.905 

112.977 

51011 

FARMV I LLE, VIRGINIA 

37 

18 

52. 043 

281 

33 

38.270 

102.434 

51013 

clearwater-st.p.fl 

27 

55 

13.006 

277 

18 

24. 529 

-32.083 

53 017 

MTFFLINVI LLE, PA 

41 

00 

57.788 

283 

39 

42,036 

256.467 

51019 

HUDSON, N. Y. 

42 

14 

37.676 

2B6 

13 

24 , 447 

12 , 477 

53 020 

AL80RG{ GSC 3 . VER. 

44 

54 

29:267 

286 

42 

30.669 

29.527 

51021 

OR L fans, NASS 

41 

51 

39.23^ 

290 

02 

59.203 

-19.649 

51022 

FAIRFIELD, ME 

44 

35 

59.522 

290 

24 

46.393 

7.587 

51023 

BOUCHARD RM2, MAINE 

47 

11 

53.91? 

291 

26 

48.671 

326.393 

51024 

FREEPORT, TX . 

29 

02 

33 - 525 

2 64. 

39 

49.664 

-31.911 

51026 

CLARKSVILLE, TEXAS 

33 

38 

2-2.877 

264 

58 

55.424 

117.625 

51027 

SPRINGOilLE, ARK. 

36 

10 

23.417 

265 

52 

39.25? 

374.242 

51028 

THftYFR, MISSOURI 

36 

34 

37.098 

268. 

22 

23.142 

259.787 

53 029 

PLATTE CITY. MO 

39 

16 

50.854 

265 

13 

49.562 

261.807 

51033 

CLAY CITY, ILL 

38 

38 

1^.569 

27] 

39 

00.669 

105.314 

51032 

EL OARA, ILL 

39 

37 

27.900 

268 

58 

33,650 

201.181 

51056 

KEARNS, UTAH 

<^0 

38 

36. 160 

248 

01 

42.569 

1369,895 

51057 

DRY 1965. NEVADA 

40 

23 

43 .717 

244 

47 

30.933 

1828.559 

51058 

DIATOM 3958, NEVADA 

39 

49 

37.575 

241 

00 

52.340 

1275.178 

51095 

AGAMENTICUS. MAINE 

43 

13 

24.265 

289 

38 

28.169 

175,472 


a = 6378145m l/f= 298.25 
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Table 3 

Coordinates in the MR-72 System (Precise Traverse) 


0 X b(i^) 


] nnoo 

CHEYFWWI' . WY. 

Af 

08' 

00 '' 23? 

255° 

07' 

57 '.' 34 1 

1889.5 

10003 

grffmwillp.imtsb. 

33 

2« 

42.387 

268 

59 

51.400 

44 „ 2 

10006 

RANCH, K/Ov. 

39 

] 3 

26.686 

261 

2.7 

29.494 

599.4 

100 IB 

jriMESTn^lN, I'EX 

30 

26 

48.204 

262 

01 

17.452 

327. 1 

10014 

FR ANKTOI'U IWO 

40 

] 4 

07.031 

274 

10 

27,148 

259.0 

100?O 

hARYWJLI.F, IwO 

3 8 

3‘^ 

20,856 

274 

21 

07. 740 

2 1 P . 0 

10021 

Cash, KY 

37 

3 3 

06. 81M 

27? 

55 

1 0 . 34 ] 

265.6 

100'2? 

IHKA, MISS 

34 

47 

15.482 

271 

45 

30,1 84 

250, 4 

100? 3 

WFRSTFR . MISS 

3 3 

33 

54 . 698 

270 

50 

04,^08 

141.3 

100?« 

FCl.TiM AFB,FL 

30 

34 

04, 2^2 

273 

47 

0] .109 

40,9 

10031 

CC’LOSTnNF, CAL I F . 

35 

25 

39.774 

243 

06 

40.716 

994,0 

10032 

FOWARS AFB.CA 

34 

57 

50.736 

242 

05 

10.754 

758.4 

1005B 

PILLAR POINT, CA 

37 


53.68 8 

237 

30 

09.319 

2.1 .0 

101 1 6 

FOWARDS AFB.CA 

34 

49 

59, 147 

242 

22 

1 3. 248 

9 39. 7 

10117 

F no ARDS AFH.C-A 


4 B 

48.361 

?4 1 

52 

23.015 

67 B. 3 

20001 

BFLTSVILLF, RiO 

39 

01 

39, 74<V 

283 

1 0 

26.713 

40 . 1 

20003 

HR ICHTHOnn.CAl 


2? 

54, 51 9 

242 

19 

09. 1 

2258.8 

2001.6 

wnniiLiNF, GA 

30 

56 

55 . 045 

27^ 

19 

07,824 

n.6 

20016 

CnLUiviHIA.MISS 

31 

12 

44. 48 1 

270 

16 

27.995 

113.6 

30025 

BinOMF TFLD, (IH 10 

40 

05 

11.744 

278 

1 5 

39.648 

360.6 

30026 

coLUNPus.nHin 

40 

00 

27.772 

276 

57 

30. 3 30 

240.2 

30027 

GR 6 FNV T LLE.OHin 

40 

09 

51,434 

275 

23 

26.755 

313.5 

3002B 

MFTAMORA, ILL 

40 

4q 

20.424 

270 

42 

40.57] 

249.0 

30029 

MOSFS LAKE, wash. 

47 

11 

07.389 

240 

39 

48.048 

355.0 

30098 

□RLAWD.CAL IFORMIA 

39 

4 A 

44 . 800 

?37 

50 

57.547 

42.6 

30099 

ChamCF 1967, MONTANA 

47 

47 

08.008 

251 

22 

07.896 

1011.9 

51008 

BOLIVIA, NORTH CAR. 

34 

02 

10.438 

281 

50 

39.492 

. 8 

51012 

BnNIFAY,FLnRIOA 

30 

39 

05.180 

274 

11 

38.234 

34.70 

5101^ 

VALKAR I A, FLORIDA 

27 

57 

25.291 

279 

26 

32.023 

1 3.258 

51015 

HI ALFAH, ELORFOA 

25 

53 

24. 936 

279 

41 

34.345 

16,3 

51025 

newton, TEXAS 

30 

54 

24. 08 0 

266 

23 

57.134 

87. B 

51030 

kimgfisher.oklah 

35 

47 

01.278 

262 

01 

12.942 

363.0 

51033 

WOODBINE. IOWA 

41 

42 

12.708 

264 

21 

21 . 168 

425,5 

51 039 

IRAAN, TX 

30 

52 

] 4.847 

258 

04 

00.426 

899.5 

51041 

ARTHUR , NEB 

41 

38 

27.013 

258 

24 

03.844 

1184 .0 

51043 

LOVELL, WY 

44 

48 

02.028 

251 

39 

19.733 

1220.35 

51044 

CRESTON.HY 

41 

36 

56. 136 

252 

12 

57.024 

2236.5 

51048 

ALBIIOUEROUE, NN 

34 

56 

43.351 

253 

32 

25.929 

1830,10 

51052 

P ICACHO, AZ 

32 

43 

24.722 

248 

29 

08.032 

472.8 

51066 

TERRE.BOMNE.OREGnM 

44 

23 

32. 167 

238 

42 

17.070 

874.9 

51067 

MINERAL WELLS.TEX 

32 

57 

44.602 

26] 

54 

36.968 

357,6 

51068 

YULEE, FLORIDA 

30 

41 

45.626 

278 

15 

59.246 

21.4 

51069 

ASHEPOO, S.G. 

32 

45 

3] . 19 4 

279 

26 

36.789 

I .9 

51070 

wrightwdod,calif. 

34 

22 

44, 426 

242 

19 

08.972 

2172.77 

5107A 

midland, DR 

42 

07 

22.217 

23« 

10 

26.461 

] 229.88 

51089 

MAYHOnO 1971, CAL IF 

38 

08 

32.353 

238 

16 

37.829 

21 . 10 

51103 

DONA ANA C0UN1Y.NM 

32 

04 

19.097 

253 

31 

06.299 

1268.1 
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Table 3 (Continued) 

Coordinates in the MR -72 System (Precise Traverse) 


hi 121 

nPFLflll.SAS. 1. 

A 

30 

3 7 

bll?3 

F(!ftT DAVIS, 

TX 

30 

40 

bl] 24 

FOEM, TX 


31 

02 

bll?.5 

MlLTnM,FL 


30 

36 

51126 

DON-TICFLin, 

FL 

30 

31 

51127 

K'OSFS LAKE, 

HA 

4? 

1 1 

52001 

bfltsvillf. 

mar Yl. 

39 

01 

62002 

ftFLTSVlLLF, 

1'' A R Y L 

3 9 

ni 

52063 

Ons PALI1S,C 

A 

36 

54 

53002 

BELTSVILLE, 

DO 

39 

01 

53063 

nos pAi.ns,c 

A 

36 

54 

60001 

bfltsvillf. 

mO 

39 

01 

9000ft 

???, Aft IZOA'A 


32 

39 


a = 6378206.4 


54, 595 

267 

50 

03,628 

22.6 

15,872 

255 

58 

38,693 

2066.7 

06.711 

259 

55 

09.670 

716.0 

06.465 

273 

02 

00.036 

47,6 

36.314 

276 

1 2 

3] .932 

55.8 

07. 389 

240 

39 

48,048 

355,0 

39.749 

283 

] 0 

26.71 3 

40 . 1 

3'^. 261 

2 «? 

10 

26 .^^99 

46 , 1 

51.1 84 

239 

26 

48.61 1 

V . 68 

39 . 26 ] 

283 

1 0 

26,899 

40 . 1 

53.183 

239 

26 

48.611 

9.68 

39.749 

283 

10 

26 . 71 ? 

4f,. 1 

10.738 

245 

24 

32.567 

61,3 

1 /f = 

294. 

CO 

00 
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Near simultaneous Geoceiver tracking data at station Nos. 51072 and 
51024 for satellite 3140 (60), over the period January 25 to February 8, 1974, 
has been received from the National Geodetic Survey along with the Precise 
Ephemeris for the satellite (state vectors at.minute intervals). 

Short Arc Geodetic Adjustment (SAGA) Computer Program had already 
been acquired earlier from DBA Systems, Inc. 

It is intended to investigate the possibility of obtaining station coordi- 
nates of geodetic accuracy with the help of on-board (and not precise) ephemeris. 
Operationally, there is generally a delay of several weeks in obtaining precise 
ephemeris, while suitable equipment (e.g. , JMR-1 equipment) can acquire the 
on-board ephemeris along with Doppler data. Therefore, if a satisfactory 
procedure for obtaining station coordinates with on-board ephemeris is formu- 
lated, it may be possible to provide a coordinated control of geodetic accuracy 
more rapidly than at pres ent . 

The investigation is proposed to be carried out in the following steps: 

(a) i) Assessment of the accuracy of station coordinates obtained with 

the available observational data and the precise ephemeris . 
ii) Theoretical formulation which may help obtain results of com- 
- parable accuracy with on-board ephemeris . 

(b) Recovery of station coordinates with the observational data and 
the on-board ephemeris, using the formulation indicated at (a) ii) 
above and comparing the results with results at (a) i). 

The data is now available on tape E 13484 (SLOT NO F 112) in four files 
as follows: 

File # 1 EPHEMERIS for Satellite ■— Station 51072 

File# 2 GEOCEIVER DATA — Station 51072 

File # '3 EPHEMERIS for Satellite — Station 51024 

File# 4 GEOCEIVER DATA — Station 51024. 
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Printout of the whole data is also available. As a result of the scrutiny 
of data, it is found that data for 30 passes is available. Preliminary runs with 
SAGA have revealed some problems in the adoption of the program. These 
problems are being looked into. 
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3. ACTIVITIES RELATED TO EOPAP 
(Grant No. NGR 36-008-204) 

3. 1 Rotation of the Earth 

The effect of crustal motions in the rotation vector of the earth may he 
studied as explained in the Fourth Semiannual Status Report , by solving the 
general Lagrange-Louiville equations under the assumption of quasi-rigiditj’’. 

The main problem then, is the obtention of [ Al]p , tliat is, the contri- 
bution to [ lol (initial va^ue of the earth's inertia" tensor) due to plate mass 
displacements . 

This value was expressed by 

[Al]p = S[Al]p (1) 

1=1 1 

where n is the number of tectonic plates constituting the earth's crust. 

The tectonic model described in [Solomon, et al. , 1975] is used in this 
work. It consists of eleven plates, of which the relative angular velocities with 
respect to a reference plate (Pacific) are given. The absolute angular velocity 
of the reference plate with respect to the underlying mantle is also given. 
Therefore, the absolute angular velocity of any plate Pi may be computed 
using the following column matrix notation 

{wjpj = {wjp^ + [wJrp (2) 

where 

{(Or }p^ = relative angular velocity vector of any plate Pi 

with respect to an arbitrary reference plate (Pacific); 

{tOa }rp = absolute angular velocity vector of the reference plate 
(Pacific) with respect to the underlying mantle; 

{tOa }pj = absolute angular velocity vector of any plate Pi with 
respect to the underlying mantle. 
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Several velocity models for [wa }rp are described in [Solomon, et al. , 1975] . 
In the present study their model B4 (continents have 3 times more drag than 
oceans) was considered, primarily because of the fact that their final values 
agree closely to the recent ones given by [Kaula, 1974] using a completely 
different approach. 

Once the absolute angular velocities [Wa }p^ for each plate are known, 
the differential changes 69 and 6X in colatitude and longitude for each 1° x 1° 
block may be computed. 

Following the theory in [Soler, 1976], this can be expressed in the case 
of a spherical curvilinear system as; 


/ 

69 


[r sin 6 cosX 

6X 

= 

H'^E,[^]p^ |r sin6 sinX 

6r 


( r cos 0 
\ 


where 

H = "matric matrix” of ihe transformation between 

spherical and Cartesian coordinates; 

R = rotation matrix of the transfoirmation between 

the geocentric and moving frames; 

[^]p^ = skew -symmetric matrix of the absolute angular 

velocity vector. 

Observe that {6co}p^ = 

After the proper values are substituted in (3), it is seen that one immediately 
obtains 


66 ~ 6 Wi sinX - 6 t 02 cosX (4a) 

6 X = 6d0i cosX cot 6 + 6 c 02 sinX cot0 - 6t0s (4b) 

where the values of {6co } correspond to the absolute angular velocity compo- 
nents for the specific plate containing the block ( 6, X ). 

Finally the elements in the summation of the right-hand side of equation 
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( 1 ) will be given by 


= S^( 60 jAlgk + aX^[Al^]k) 


(5) 


where m is the total number of 1° X 1* blocks on the plate Pj and 


[A 10]!. 


s • differential changes in the earth’s tensor of inertia 
due to a differential motion 66 k of block k; 


[Al^lk = differential changes in the earth’s tensor of inertia 
due to a differential motion 6 Xk of the block k. 


The computation of [Al 0 lk and [Al^^^lk involves the integration over 
the volume of every block k, taking into consideration the theory of isostasy 
according to Heiskanen. 
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3,2 Close Grid Geodyna-mic Satellite (CLOGEOS) System 

The preliminary investigation with a grid network of 9 stations at 5 
minute spacing in both latitude and longitude in the San Andreas fault zone 
for the Close Grid Geodynamic Satellite (CLOGEOS) system was completed 
during the period under report. The results were published under Reports 
of the Department of Geodetic Science No. 230, OSU. 

A more detailed simulated study with 75 stations along the three 
major faults (viz. , San Andreas, Hayward and Calaveras) in central Cali- 
fornia has been started. In addition to the experiments conducted in the 
preliminary investigation, the present effort will include the following; 

(L) Simulation of weather effect and actual crustal motion; 

(ii) Study of the effect of observation pattern or station grouping in 
geometric and short arc modes; 

(iii) Study of the effect of significant digits in computing a solution in 
near critical station configuration; 

(iv) Geodetic monitoring of cmstal motion. 

The results are expected to be included in the next status report. 
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3. 3 Modeling of VLBI Observations 


Let two stations denoted i and j be engaged in a VLBI observation of a 
radio source p at some epoch t^. The observed time delay t of arrival 
of a certain wavefront at the two stations provides a distance d^pk = c ‘ t 
which is the projection of the station-to- station vector By on the instantaneous 
unit vector of the radio source direction ep*(tk). Both Bij and e^ (t}.) refer 
to a reference frame fixed to a network of stations in a certain way. Obviously 

^IJpk “ (^k ) • 

If Cp is the unit vector in the direction of tlie radio source with respect to an 
inertial (source- fixed) frame, then one has 


8p*(tk) = M(t^)ep. 

The transformation matrix M can be parameterized in terms of Eulerian 
angles as 

M = R3{'^) Ri(6) R3(!/>) • 

If a deterministic model for ttie time variation of 0, ij) were available 
in the form 

d^ E 

= f(E,t) E = [<p, 0. 

then solving this equation, one should have E(t]j) = F[E(to ), E (to) ] , or 

M(ti_) = M(<po, 00, ^0, 00 ; ^0, tfc) 

with a total number of six parameters . However, such an approach is not 
possible because of the uncertainties surrounding current laaowledge of the 
earth’s rotation. Alternatively, one may set 

M = Rs{(P'’+5(D) Ri(0’’ +6 0) R3(!/)° + 6j/j) 

where 6°, fj)° are approximate values provided from the analysis of 
classical astronomical observations and 6<p , 66, 6 are small corrections 
which are assumed to be constant over a short interval of time (e.g., one day) . 
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The niunber of parameters in M is now only 3 , as compared to the previous 6. 
However, the direction and magnitude of the instantaneous vector of rotation 
00 is one of the objectives of this analysis, and since 

CO = co((??,e,0,(p, 0,0) 

one has to interpolate between values of <p = + 6 <p , 6=G° + 6 0, 0=0° + 60 

to obtain 0 , 6 , 0 . This additional approximation destroj^s the optimality in 
estimates of to . 

For this reason, and despite the fact that 3 parameters are sufficient 
for the description of the relative rotation of the two reference frames, a 
6 parameter model is used which explicitly contains the co vector, 

M(t,) = Rier?)R 2 (-^) R3[n(tk-to) + 0o] R 2 (E)Ri(H) 

where Q is the angular velocity of the earth’s rotation and the geometric 

meaning of the rest of the parameters is depicted in Kgures 1 , 2 and 3 , 

Also, ^ and tj are the usual coordinates of polar motion; H and H are 

two similar parameters describing precession -nutation and 0o + f2{t^ - to) 

is an analog of GAST (Greenwich Apparent Sidereal Time). A priori 

— ^ 

approximate knowledge of the co direction with respect to both inertial and 
network frames, can lead to selection of coordinate systems so that 
E, H, ^ , 17 are small quantities with zero approximate values. Then 
one can write to ahi^ degree of approximation 

/I 0 M no -E 

0 > H 

Vl 1) I / -H I 

Setting Bij = (Xj - Xj, y^ - yj, Zj - Zj)% tj. = tj. - to and Cp = (cos6p 
coSQTp, cos6p sinttp, sin6p), the model becomes 

t 0 |\ f\ 0 /cOS^P a>SKp\ 

^ J -n i H sinotp 
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The partials with respect to parameters evaluated at approximate values 

= become after setting “ ©o + Tk > Xkp~ ~ cep» 

Xij = Xj - Xi, Yii = Yi- Yi and Zjj = Zj - : 


M. = 

31 




5 d _ 
9 n 

5 d ^ 
3H 

3 d ^ 
30 , 

M - 
3 a 

3 d ^ 


- 2 SJJ cos^p - -^y sl-nlTp 

(-«ij cosy^ + stny^) sinS’p + Sy cos«Xp 

( jfig- sinyp + cosy, ) sm^'r “ cos^p Siiaecp 

-(“ij ®'’'aKp+ 

•p 3 d. 

“ 300 

_ = OOS^P ^®^XkP 

3 St ^ 


3 d _ „ _ _ cos^p 

9 ^i 

^4 _ _ 34 _ sm S'p 

|i = (-X;- sinx^p + «sx,p) cosf, 

|i „ (-X;. «$x,, + s'mxu,) slnSp + a cosSp 

3 op 
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3. 31 Definition of Coordinate Systems, Inner Constraints 


Linearization of the observations leads to the equations 
L = AX+ V 

where L contains differences between observed and computed distances; 

A contains the partials already derived; X contains corrections to parameters 
and V is the vector of random residuals. Setting the weight matrix to 
P = I for simplicity, the normal equations become 

(A^A)X = A^L or NX = U. 

But because of a lack of system definition, one has rank (A) = ?7 - 9, where 
r) is the number of parameters and consequently N is singular. The rank 
deficiency 9 of A is due to the 9 degrees of freedom of the undefined coor- 
dinate systems: 3 for the orientation of the inertial system, 3 for the 

orientation and 3 for the origin position of the earth-fixed system. 

Coordinate systems can be defined by means of 9 minimal constraints 
on the station coordinates (dxj , dy^ , dzj ) and the radio sources coordinates 
(dap, d6 p) . The number of constraints can be reduced if use is made of 
naturally defined observable directions — in the present case — the direction 
of the earth's vector of rotation co . One may then assume that the Z axis 
of the inertial system coincides in direction with O), (but not with the.z axis 
of the earth-fixed system), and tlius avoid or limit the variation of station 
coordinates with time. This leads to elimination of the columns of the design 
matrix A corresponding to the parameters E, H, so that onl 5 '' 7 constraints 
are now necessary. One has rank (A) == rank(N) = rj - 7. 

Among the solutions to the normal equation, the unique one given by 
X = N U has the properties X X = min and trace N = min. hi view of 
the interpretation of N as the variance -covariance matrix of parameters, 
the second property malces the solution optimal. To avoid the use of a 
pseudoinverse computation algorithm, one may construct a set of inner constraints 
giving the same solution as the pseudoinverse. The inner constraints can be con- 
structed with the help of the geometry of tlie operator represented by the matrix N. 
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The domain JD of this operator is the space of all n x 1 column vectors and 
can be turned into a complete inner product space introducing the inner product 
<f»g> = g^f. Then JD has the ox'thogonal decomposition 

J) = X© 


whereX is the null space of the operator, i.e. , the set 


{y s J); Ny = 0 ] . 

It is well known that X = N U e X'^i and therefore the condition X J- X 
uniquely defines X, If [Ci 3 is a basis for X" (i 2, • * • 7, sinceX^ is of 
dimension 7), then the condition written 


Xj-Cj 

In matrix notation: 


<X, Ci>=C,X=0 i = i.2,..,7 


X = (c, c^...c,) X = c'^x = 0 , 



The problem thus reduces to finding a basis inX? i- e. » in finding 7 linearly 
independent n X 1 vectors Ci satisfying 


NCj =' 0, or since N = A^A, ACj = 0. 


If Aj is the jth row of A, one must find 7 linearly independent solutions 
y = Cj, 1 == 1, 2, • • ■ 7 to the set of equations 

A^ y = 0 j = 1, 2, • • ■ s (s = number of observations) . 

Setting y' = 3 


where m = number of stations and k = number of radio sources . 

For the row of A corresponding to an observation dij-pn (stations i and j observing 
radio source p at epoch t^), one obtains 




^ + g -M. 0 

aotp as*. 
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Making use of the analytical expressions for the partials and after a con- 
siderable computational effort, one arrives at a set of solution vectors 
{Cj which gives rise to the inner constraints C ^ X = 0 with the matrix 
being 

dQ^ dn dx, 

ObO 0 ^ 00 I OQ...I 0 0 

000 0 0 10 0 tO-'-O > 0 

006 0 OOl 0 0 j 

1 0 0 0 0 -a, 0 -Xg . . 0 

0 I 0 0 5, 0 2.g 0 . • 2^ 0 

0 0 -1 0 X, 0 “1^5, 0 

0 0 0 0 0 C 0 

I 

■det, d§^, d«^ 

I 
I 

JO 0 0 0 b 

I O' b O 0 

Jo 0 0 

I 

I 0 0 0 0 . • > ' • 0 

!' 0 0 0 6 . • ' ’ • Q 

I 

'0 0 0 0 ■ • ' ' ‘ 0 

1 I 0 t 0 . . . • . t 


More explicitly, there are 3 sets of inner constraints. The first defines 
tlie origin of the earth-fixed system 

dsfi 

dt^i = 0 . 

dZi 
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The second set defines the orientation of the earth- fixed system 




Finally, the third set defines the direction of the X axis of the inertial frame 

k HI 


i 


P=i 


t=i 


3. 32 The Hole of Minimal Constraints in System Definition for 
a Non-rigid Network of Stations 

The solution for station coordinates using minimal constraints depends 
on die set of used approximate values. The motions of stations with respect 
to each other can be modeled as follows: Consider a sequence of time inter- 
vals Ati == [ti_i, tj] i = 1, 2, • ■ • , and assume that within each interval 
that station coordinates remain unchanged. At the epochs t^ the coordinate 
values "jump" to a new set of values, i.e., the variation of coordinates with 
time is modeled by simple "step functions. " Now observations within each 
interval Atj can be treated separately in a sequential fashion. If for the 
Ati interval one uses as approximate values of the coordinates, their esti- 
mates provided from the analysis of observations in the interval Ati_x, then 
the minimal constraints provide a means of system definition for a deformable 
network of points . 

The coordinate axes thus defined are "geographical axes, " i. e., they 
are presex’ibed to the station in a specified way. From this point of view 
inner constraints may be inappropriate because they involve the direction 
of the rotation axis of the eartli. Some more appropriate and intuitively 
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appealing sets of minimal constraints are discussed next. 

For the position of the origin of the system,' the constraints 



already contained in the inner constraints assures tloat the coordinates of 
the center of mass of the stations (considered of unit mass) will remain 
the same. For the definition of the direction of the system axes, two choices 
appear to be dynamically meaningful . 

One is to consider the axes of zero relative angular momentum for the 
stations unit mass points in a fashion similar to the definition of the Tisserand 
axes for the whole earth. If rj denotes position vector of the ith station, then 
the total relative angular 'momentum vector of the network is 

=0 or n ° 

V di i 

In matrix form one has 


■m 

0 -2i 



z 

Zi 0 'Jfi 


H'i 

i~\ 



dZi 


A second choice is the principal axes of the network, i.e. , the choice of 
axes that malces products of inertia of the stations vanish. To retain the 
zero (or some other constant) value of these products one must set their 
variation equal to zero . The products of inertia are 


j t i ’ *7- 

dr> = r ( S. + = 0 , 


i 


or, in matrix form 



0 

3 .; 


S; 







dZ{ 






= 0 
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Rapid City Journal 




Land use 


By DICK REBBECK 
Journal Outdoor Editor 
STURGIS — By no means unique 
in its problems of subdivision 
sprawl, Meade County is welt into a 
comprehensive land use planning 
demonstration that promises lo 
develop ideas other counties might 
■well apply to their situations. 

'Describing the problem for the. 
project, Arnold Bateman, rural de- 
velopment speclal.lsc for 
South Dakota State University, . 
has" specified "uncoordinated 
growih occuring adjacent to the 1- 
90 corridor and contiguous to the 
city of Sturgis.” 

Kirk Carlsten, Meade Coumy 
planning and zoning coordinator, 
points out that most of (he county is 
agricultural, with little prospect of ' 
residential or Industrial impact:’ 

‘ ’99 per.cent of the problem is along , 
!*90, from Ihe (Black Hills Nation* 
al) Cemetery to Black Hawk." ■ 

* ‘'Much of the development 
(there) is within steeply sloping, 
forested areas, which increases 
.er(^sion potentials; creating the 
possibility of stream slltation and 
pollution or water source contain* 
mation," Bateman reports. 

. “Such development presents 
problems of sewage disposal and 
|adequaie ' watdT^ppl/ ' foT "all • 
\dcmands, Including firefighting re* ) 
‘qulrements, congestion which may > 
-.jcxcecd. (ha-carrying capacity of tha • 
’lank, increased demands for police f 
and fire protectiQn.ond other public ^ 
services, such as adequate roads," 
he explained. 

Similar land use questions have 
aHsm^ adjacent _to_ElJsworth Air 
■“’Porce^Base'andrnenr'B 
State Park. i ’ ’• • - ' 

j , ■ ?•' _ 

At Bear- Butte,. Carlsten ’notes, 
present 20«acr^ ownerships are pot 
within the scope of a newly adopted 
subdivision ordinance. 'Tf some- 
one wants to subdivide a 20-acre 
site. Il might come under- it, he’ 
says. 

• Several subdivisions have been 
'well laid out, from a land use p)ann>' 
ing aspect. Evidence of this, cited 
by Carlsten, includes community 
water systems, which avoid san- 
itation problems from septic tanks, 
provisions for safe waste disposal, 
provisions for drainage .andJ^lODd 
protection, access for emergency 
vehlcJes, water supply* for fire- , 
Fighting and grouping of likc-vdlue 
houses. 

"In residential zoning, wc-should 
aim for more than one type of res- 
idential zone, and in each area . 
define the size of lot, ^uare foot- 
age of buildings and some other 
covenants," Carlsten comments, 
‘Tf people put good money into a .. 
place ... the value of their proper- 
ty goes down if someone else puls a 
lower value house in with the 
* higher ones." 

All gradations of building can be , 
'provided for through zoning vdth- • 
out jeopardizing Investment, he 
believes. ' 

Mobile, homes and modular 
homes present a special case in this 
context. Phil Cervany,^ physical 
planner with the sixth District 
Council of Local Governments, 
’says mobile, homes historically 
have been confined to commercial 
zones because lots were rented for 
profit. But it’s actually a residen- 
tial type of use. and one that's 





^Spectacular views make Sturgis area subdivisions popular (Journal photos) . ' 


growing as economic factors die- the whole planning- operation is to 
tale greater use of, this type of provide in every respect a belter 
•housing. Zoning should provid^for environment for the people,’’ Bate- 
It as a type of residential zopc7 he man reasoned, "the citizens of the 
indicates, county must be given the opportun- 

j liy to take part in the planning 

Observing one mobile home area process." 
where residents buy their one-acre Bateman further emphasizes, that 
lots, Cerveny suggests ownership goals' and objectives must be the . 
motivates people to take better primary responsibility of citizens ' 
care of their properly. and tlieir local officials "and not 

the responsibility of the technical 
^ While acknowledging good land planner." The planner implements 
'Use judgments in some develop- within guidelines laid down by the 

ments, planning, specialists and people. ^ . ’ 

county of Flcials have detected prob-. 

lems along the 1-00 strip, not the^ In Meade County. the_pdncipal 
- leastof which was tbe joss of prime . . -goals -and objectives.have been'd^ 
agricultural lank.’ terfnlned to be; 













im 


Carlsten finds rernote sensing map valuable tool 


• ’ “The inevitable deterioration of 
the . . , natural environment, quali- 
. ' ty of living and the increasing de- 

■ mand'jor pbblic services became,, 
apparent and the Meade County 
commissioners began looking for 

• solutions,” Bateman reported in an 
'] Agricultural Experiment Station 
article soon to be published. 

Out of such concern came the 
.countywide comprehensive land 

■ use planning pilot project, one 

grounded in local public involve- 
ment. ^ 

"Since the general objective of 


ling- operation is to in response to 1974 state leglsla- 
y respect a belter tion, the county se^ out to 'set up a 
■ the people,’’ Bate- comprehenslyq land use plan by the 
"the citizens of the July i,. 10':^$, deadline. Elements in 
given the opportun- this plan are to include a zoning or- 
rt in the planning dinance, a planning ordinance and a 
■ land use map. 

her emphasizes, that . 

ctives must b4 the A county subdivision ordinance 
isiblUty of citizens ' 'tiht went into effect Jan, 1 essen- 
otficioJs "and not tlally provides for "harmonious de- 
ty of the technical velopment” through coordination 
dunner implements different developments so that, 

■s laid down by tbe hi Cerveny's words, an undue ft- 
, nanclal burden won’t fall hack on 
" ' ‘ . -the county” for such public ser* 

ounty. the principal vices as road mumtenance schools ■ 
ctives.have been-ap- -“h" fire ‘‘'><1 Puli':® Pratoetton. 

; Weather and high Interest rales 

. have-* "pretty well brought con- 
' struclion starts to zero,” sb-tarlj.' 

Bten doesn’t have a measure of how. 
— ^the subdlvislon- ordinance will func-' 
tion. "We won’t have a good eval- 
nation until spring and summer." 

^ Meade County is now working on 
‘ Its eoumy, -zoning- ordlhaaca. ' Carl— 
believes Meade County Is as 
along toward the 1976 deadline 
as any West River county. The zon- 
ing ordinance Is expected to be 
ready for the county commlsslon- 

monihs,” "I don't see building 
i ; 1^. codes for some time to come. In tbe 
county," he adds. "That will take 
quite an enforcing arm. The cltyils" 
county is not." 

* ' V, ^ Work is also proceeding on a soils 

\ f - ^ maip using aerial photographs ob- r 
" ' f tained in cooperation with the Re- \ 

'v mote Sensing Institute at I 
- ^ '• / Brookings. This will be especially A 

•betiei^agrl-*' 

I cultural jand, something people of 
r' county have indicated they 

've zone an area strictly ag-- , 
rlcultural, we’l) protect it totally,” • 
Carlsten says. "If it gets zoned res- , 

■ ■ idential and it’s in agricultural use, 1 
restricted for agricultural 
agricultural). But if it Is j 
' subdivided, it will be considered h 
residential in all senses.” 

^ landowner still has the right to 
’ subdivide, m a residential zone, and 
. ^ ' the county doesn’t have to go 

he land base sup- through rezoning II from agrlcuJ- 
Iture, forests, natur- tural use. At' the same time, 
and .minerals; with though, it is hoped that this\will en- 






® Protect the land base sup- 
porting agriculture, forests, natur- 
al resources, and .minerals; with 


orderly land use change to meet » courage retention of someagncul- 


people's needs and protect key sites 
for future development. 

« Preserve a high quality envi- 
ronment through pollution' abate- 




protection, maintaining recreation 
the and other special-value areas and 
and avoiding particular hazards, siich j 
one as 'flood plains. . 

a Encourag? development of^* 
pleasant, efficient and safe cotti 4 
3 cf mumties. ^ 

spiiMis)iiuGi]^naW:OE^^ 

A M.'fa 


(UTdl and open space land within 
the I*9D strip. 

The land use map will also heli>li-, 
identify flood bazard.,..hi,gff*'waterl 


ment, soil erosion control, wjldljfe j table, .soil- slippage and other prob- 


lem areas lo be accounted for in | 
any future development. 

‘Tm. sure we'll have the 1^4 i 
with which we can do a better job,” j 
Carlsten summarizes, "If we don’t 
do a better job, It’ll be our own 
fault.” 

pi 


"GEOS-I SECOR Observations in the Pacific (Solution SP-7)” 

National Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 7-10, 1970. 

"Investigations of Critical Configurations for Fundamental Range Networks" 
Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C., 

April 15-17, 1971. 

"Gravity Field Refinement by Satellite to Satellite Doppler Tracking" 

Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C. , 

April 15-17, 1971. 

"GEOS-I SECOR Observations in the Pacific (Solution SP-7)" 

Symposium on the .Use of Artificial Satellites for Geodesy, Washington, D, C. , 

April 15-17, 1971. 

"Separating the Secular Motion of the Pole from Continental Drift - Where and 
What to Observe ? " 

lAU Symposium No, 48, "Rotation of the Earth, " Morioka, Japan, May 9-15, 1971, 

"Geodetic Satellite Obsei'vations in North America (Solution NA-8)" 

Annual Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 6-9, 1971. 

"Scaling the SAO-69 Geometric Solution with C-Band Radar Data (Solution SC 11)" 
Annual Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 6-9, 1971. 

"The Impact of Computers on Surveying and Mapping" 

Annual Meeting of the Permanent Committee, International Federation of Surveyors, 
Tel Aviv, Israel, May 1972. 

"Investigations on a Possible Improvement of Terrestrial Triangulation by Means 
of Super-Control Points" 

lAG International Symposium - Satellite and Terrestrial Triangulation, 

Graz, Austria, June, 1972, 

"Free Adjustment of a Geometric Global Satellite Network (Solution MPS7)" 
lAG International Symposium - Satellite and Terrestrial Triangulation, 

Graz, Austria, June, 1972. 

"Conjugate Gradient Method (Cg-Method) for Geodetic Adjustments" 

Annual Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 3-6, 1972. 



"Preliminary Results of the Global Satellite Triangulation Related to the NGSP” 
Journees Luxembourgeoises de Geodynamique, Luxembourg, February 19-21, 1973 

"Present Status of Global Geometric Satellite Triangulation and Trilateration" 

54th Annual Spring Meeting of the American Geophysical Union, Washington, D. C. 
April 16-20, 1973. 

"Free Geometric Adjustment of the OSU/NGSP Global Network (Solution WN4)" 
First International ^mposium on the Use of Artificial Satellites for Geodesy 
and Geodynamics, Athens, Greece, May 14-21, 1973. 

"Earth Parameters from Global Satellite Triangulation and Trilateration" 
International ^mposium on Earth's Gravitational Field and Secular Variations 
in Position, Sydney, Australia, November 26-30, 1973. 

"Review of Problems Associated with Geodetic Datums" 

International ^mposium on Problems related to the Redefinition of North 
American Geodetic Networks, Fredericton, N.B., Canada, May 20-25, 1974. 

"Marine Geodesy - Problem Areas and Solution Concepts" 

International ^mposium on Application of Marine Geodesy, Battelle Auditorium, 
Columbus, Ohio, June 3-5, 1974. 

"Station Coordinates and Geodetic Datum Positions froni the National Geodetic 
Satellite ProgTam" 

First Pan American Congress and the 

Third National Congress of Photogrammetry, Photointerpretation and Geodesy, 
Mexico City, Mexico, July 7-12, 1974. 

"Review of Classical Methods for the Determination of Geodetic Datums" 
International Colloquium on Reference Coordinate Systems for Earth Dynamics. 
(lAU Colloquium No. 26) 

Torun, Poland, August 26-31, 1974. 

"Global Satellite Triangulation and Trilateration Results" 

Mtercosmos Symposium on Results of Satellite Observations 
Budapest, Hungary, October 21-24, 1974. 

"Crustal Motion Monitoring with the Proposed Close Grid Geodynamics Satellite 
Measurement System" 
lUGGXVI General Assembly 
Grenoble, France, August 16-31, 1975. 
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"Western European Satellite Triangulation (WEST) Station Coordinates in the 

OSU WN14 System" 

lUGG XVI General Assembly 

Grenoble, France, August 16-31, 1975. 

"Aspects of Positioning Using Satellite Borne Laser or EF Systems" 
Symposium on New Methods of Space Geodesy 
Leningrad, USSR, November 24-30, 1975. 
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Attachment 1 



U.S. DEPARTMENT OF C0P,/J1VJERCE 

National. Oceanic and Atmospheric Administration 

NATIONAL OCEAN SURVEY 

Rockville. Md. 20852 C13/BKM 


September 16, 1975 


Dr. Ivan I. Mueller 
'Department of Geodetic Science 
The Ohio State University 
1958 Neil Avenue 
Columbus, Ohio 43210 

Dear Dr. Mueller: 

In accordance with our telephone conversation about two months ago, 
data for 55 stations along the transcontinental traverse network 
are enclosed. These data include NA 1927 Datum positions, elevations 
and geoid heights; M-R ’72 Datum positions (see comments attached); 
and Doppler results for each station. 

The data sheets are arranged in numerical order by station number. 

A sketch of the transcontinental traverse net, attached to the data, 
shows the approximate location of each station. Please note that 
Doppler positions were redetermined at the following four stations. 

1 - 20001-52001-60001 

2 - 52002-53002 

3 - 52063-53063 

4 - 30029-51127 

The M-R ’72 Datum positions and geoid heights are considered as 
preliminary, A simultaneous adjustment of the T.T. net will be 
performed after the field surveys are completed. The expected 
completion date is July or August 1976. 

Sincerely, 

B. K, Meade 

National Geodetic Survey ’ (Retired) 

Enclosures (2) 
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Comments on M-R '72 Datum Geographic Positions 


‘he geographic positions of stations given on the Doppler 
lata sheets, identified as the M-R '72 -Datum, v/ere obtained 
Torn adjustments as follows. 

1 - Western loop adjustment of the transcontinental 

traverse. This loop involves stations 10006-10018- 
51103-30098-30099-10006. The NA 1927 Datum position 
of MEADES RANCH was used for position control. 

2 - Eastern loop adjustment of the transcontinental 

traverse- This loop involves stations 10006-10019- 
20001-51068-20016-10018-10006. Also the loop 
involving station 10003 and section from 20016 to 
junction near 10019. The NA 1927 Datum position of 
MEADES RANCH was used as position control. 

3 - Northeastern section of the western loop from junctior 

north of MEADES RANCH to 51044 to 30099. The junction 
point north of MEADES RANCH and station 30099, as 
determined in the western loop- adjustment , were used 
as position control. 

4 - Positions of stations 51014-51015; 51048; 10031; and 

10055 were, determined from spur adjustments with 
control from the main traverse loops. 

5 - Stations 10018, 51067, and 51030 are common to the 

eastern and western loops. The positions given for 
these stations are the mean values of results from 
the eastern- and western loop adjustments. 


Comments on Doppler Data 

•The original version of NOAA Form 76-178 gives the- Doppler' 
X-Y-Z coordinates and height above the elli^oid referred- 
to the tracking equipment reference point. ^The height of 
the reference point above the mark is given v/1th the data.^ 

The Doppler positions and ellipsoid heights furnished by 
DMATC, Form 115-84, have been corrected to agree with 
results computed by DMAAC, dated October 7, 1974. ~~ ^ 

All Doppler data to be used in the new adjustment of the 
NA 1927 Datum, determined from computer programs of other 
organizations, will be recomputed using the NGS program. 

B. K. Meade 

National Geodetic Survey (Retired) 

September 16, 1975 
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* Doppler positions to be reobserved during 
the summer or fall of 1975. 
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U. si Network of transcontinental traverse and Doppler 
stations at points along the traverse. Doppler positions 
have been observed at additional traverse stations along 
the interior sections of the western loop, however, 
adjusted positions based on the M-R ’72 Datum are not 
available on this date, 19 75, 
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Attachment 2b 


Western European Satellite Triangulation (WEST) 

Station Coordinates in the OSU WN14 System 

Alfred Leick, M. Arur and Ivan I. Mueller 
Department of Geodetic Science, The Ohio State University, Columbus 


1. Introduction — The Ohio State University acquired the data collected during the WtEST program 
to improve the values of some station eoordiirates on the European continent which are presently 
included in the OSU WN14 solution [3], The secondary aim of the solution was to add some new stations 
and to assess the quality of the WN14 solution with the help of the additional data available. 

The WEST optical data was available in two form’s . In the first form the data comprised tiie 
direction cosines of a sin^e fictitious image per plate and the corresponding standard deviations 
derived from polynomial fitting. The second form contained the direction cosines of seven fictitious 
images without statistics. 

Solutions with flie single image data have been completed and the results are summarized 
here in a brief form. More detail may be found in [6] . 

2. Method of Computation and Results — A modified version of OSUGOP [5] was utilized to 
obtain the normal equations .and to perform the adjustment . 

Table' 1 gives a list of all the stations for which observational data was available. The 
stations had to be renumbered in order to avoid confusion with the WN14 numbering system. 

Stations which already appeared in the original wrN14 solution were constrained in the new 
adjustment to their WK14 coordinates (see (3j) with weights compatible to their a posteriori 
variances. These stations are indicated with an asterisk in Table 1. 

For appropriate scaling the chord 6C06 TEOMSO-GOie CATANIA was constrained to 
3,545, 871.454m with a weight corresponding to' 1: 10®. 

Relative constraints were also applied to maintain the relative positions of nearby stations . 
These relative constraints are based on survey information available in [1] , [2] and [4] (Table 4.) 

Ellipsoidal height constraints were also applied after transforming the European Datum 
height information available in [2] to the WN14 system (Table 3.)’ 

The final coordinates as a result of the solution with, the single image data are presented 
in Table 1. 

The transformation parameters between this solution and other recent solutions which give 
the coordinates of common stations are presented in Table 2. 

3. Conclusions — In all cases the variances of tire ex-stations in fee \VN14 solution have 
improved by utilization of the WEST single image data. The coordinates of the stations 8705 
(BRDUX), 8712 (OPICI), 8713 (ORL4A) and 8714 (SRDIN) still erfiibit extraordinarily large standard 
deviations. This is an immediate result of increased variances in the observational data. See [2] , 
Table 4, which gives the standard errors stationwise as obtained by the smoothing pi'ocedure. 

The adjusted coordinates of 8711 (CATAN) were also expected to exhibit a large standard deviation 
according to [2], Table 4, but this station has been connected to the nearby WN14 station GOIG 
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(CATANIA) by relative constraints, and therefore does not exhibit a large variance. The large 
standard deviations at station 8722 (REICVK) are due to unfavorable geometric conditions and 
insufficient mnnbers of observations. 

The transformation parameters in Table 2 show that the origin of the present WEST 
solution coincide witli the 1W14 origin, as was intended. The parameters for the ED 50 dafui 
can be regarded as reliable, while shifts with respect to the solutions NGS, SE3 and GEilG 
should not be overemphasized due to the lack of adequate numbers of common stations. 
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Table 1. WEST (33), Coordinates in the WN14 System 
(all units in meters) 


Solution 

KGS 

(comb.) 

SE3 

WN14 

GEM6 


No. of Stations 

3 

5 

8 

CO 


AX(m) 

-11.41:1.8 

-11.9±3.5 

-0.4±1.2 

-22.5±2,3 

-96.6±3.3 

AY(m) 

-10.6±2.2 

-20.4±4.1 


-30.1±2.9 

-122.US.4 

AZ(xn) 

11.6i2.0 

9.2±3.5 

-0.6il.4 

6.1±2.5 

-126.31:3.2 


Table 2 , Transformation Parameters (WEST (33) - Other System) 
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Table 3. Height Constraints* 



Ellipsoidal 



Ellipsoidal 


Station 

Height (m) 

^(m) 

Station 

, Height (m) 

^m) 

6006 

113.2 

4.0 

8708 

186.7 

3.0 

6016 

16.3 

4.0 

8709 

949.3 

3.0 

6065 

960.1 

2.5 

8710 

80.4 

3.0 

8009 

41.1 

4.0 

8711 

15.4 

3.0 

8010 

920.6 

2.5 

8712 

393.4 

3.0 

8011 

135.0 

4.0 

8713 

194.4 

3-.0 

8019 

394.7 

4.0 

8714 

144.4 

3.0 

8030 

183.2 

2.5 

8715 

15.4 

3.0 

8701 

490.8 

3.0 

8716 

686.8 

3.0 

8702 

114.1 

3.0 

8717 

686.8 

3.0 




• 8718 • 

143.7 

3.0 

8703 

51.4 • 

3.0 

8719 

1740.9 

3.0 

8704 

164.5 

3.0 

8720 , 

40.8 

3.0 

8705 

107.3 

3.0 

8721 

300.6 

3.0 

8706 

224.2 

3.0 

8722 

. 0.6 

3.0 

8707 

81.1 

3.0 

8723 

13,5 

3.0 


*with respect to the WN14 system and a = GSVSISB. Om, b ='6,356, 769.7 m. 


Table 4. Relative Position Constraints 


From Station 

To Station 

A^(m) 



8711 

8719 

13329.91 

10072.67 

-22958. 92 

8711 

8715 

-4.04 

-8.24 

7.88 

8715 

8719 

13333.95 

10080. 89 

-22966.80 

8707 

8710- 

1.78 

-1.33 

-0.16 

8716 

8717 

-0.06 

2.70 

0.28 

8009 

8723 

3709. 06 

1053.54 

-2925.82 

8711 

6016 

-1.61 

-0.43 

2.17 

6065 

8709 

20.79 

53.37 

-25.09 
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Attachment 2c 


Crustal Motion Monitoring with the Proposed 
Close Grid Geodynamics Satellite Measurement System 

Ivan I. Mueller, B.H.W. van Gelder and M. Kumar 
Department of Geodetic Science 
The Ohio State University, Columbus 

Abstract 

The Close Grid Geodynamics Satellite Measurement System consists of a 
satellite -borne laser or microwave system and closely spaced inexpensive reflectors 
or transponders on the ground. The feasibility of the system is investigated and 
the expected accuracies, based on extensive simulations, are presented. Various 
possible applications, including crustal motion monitoring, are discussed. 

* 

1. Introduction. -- Originally, CLOGEOS was conceived as an orbiting rang- 
ing device with ground base reflectors. A grid of these reflectors (spacing 0.5 - 

50 km) which are projected to be low cost (passive, maintenance free and unattended) 
will permit the saturation of a local area to obtain data useful for geodynamic and 
• geodetic (oceans included) purposes. In this investigation a first attempt was ‘made 
to get an insight on how maximum accuracy of relative station positions can be 
achieved in a short time span (3-5 days). 

2. Instrumental Concepts . -- Measurement systems as laser radar, RF 
radar or a combination of both operating in continuous wave or pulse mode are 
able to provide ranges, range rates (Doppler) or range differences (integrated 
Doppler). 

In this study only ranges were considered with the already feasible laser pre- 
cision of 10 cm. The ranges are observed in two modes, ‘ simultaneous and non- 
simultaneous . 

Two types of vehicles carrying the transmitter have been considered: A. 

Satellites at various altitudes: 392, 657, and 1007 km. The satellite orbits (passes) 
were generated with the Goddard Trajectory Determining System (GTDS), developed 
at NASA's Goddard Space Flight Center [1>2]. B. Airplane flying at an altitude of 
9 km. 

3. Ground Stations , Two types of stations were considered 

(Fig. 1): A. Nine grid stations with a spacing of 5' were chosen in the vicinity 

of the San Andreas Fault area in California ( Ato = 9*3 km and AX =7.3 km). The 
ellipsoidal height differences between the stations were varied between 0 and 1000m. 
B. Three distant reference stations were selected outside the grid area near 
San Diego, and Quincy in California, and near Bear Lake, Utah. - 

4. Recovery of Relative Positions of Grid Stations . -- Having simultaneous 


Presented at the XVI th General Assembly of the lUGG, Grenoble, August 18- 
September 6, 1975. 
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and non-simultaneous ranges two different algorithms can be used to compute the 
relative positions of grid stations: A. Geometric adjustment which takes advantage 

of the simultaneity of the observations. The software used was the Ohio State 
University Geometric and Orbital (Adjustment) Program for Satellite Observations 
(OSUGOP) [ 3 ]. B. Short arc adjustment (dynamic mode) which does not have the- 
requirement of simultaneous observations. The' software used was the Short Arc 
Geodetic Adjustment Program (SAGA) [4,5], 

Since a range measurement system lacks any coordinate system definition, 
especially in the geometric mode, the recovery of the relative positions was express- 
ed in terms of the estimable quantities, the lengths of the chords between the grid 
stations (Fig. 1} and the angles between the chordsC 6 , 7 l. 

5. Georhetric Mode Results . -- The geometric mode leads to a very simple 
mathematical model. However, local satellite ranging networks often degenerate 
into critical configurations (see Table 1, line 1) as opposed to global satellite rang- 
ing networkstS^.To avoid these critical configurations two possibilities are mentioned: 

1. Separate stations in height either by giving the grid stations some height differ- 
ence AH (Table 1, line 2), a possibility only in the case ofi accomodating topography, 
or by including into the observation campaign the three reference stations outside 

the area (Table 1, line 4), This possibility has the stringent requirement of having 
favourable weather conditions at 4 different sites (grid area and 3 reference stations). 

2. Separate the ranging devices in height. The best (and most realistic) solution, 
to avoid the effects of critical configurations within the limited area of the grid is 
the combination of an airplane and a satellite (Table 1, line 5, iSTote that no distant 
[reference] stations are needed). The only disadvantage of geometric mode is 

the instrumental problem related to the realization of the simultaneous observations. 
These at least for the lasers may be overwhelming. 

6. Short Arc Mode Results . -- The absence of the requirement of having 
simultaneous observations and' the absence of the bothersome critical configurations 
are the main advantage of the short arc mode. However, in order to get stability 
in the solutions the 3 distant (reference) stations must be observed during each 
pass (Table 1, line 6). A pass of 4 to 10 minutes lengths for satellites at altitudes 
between 400 and 1000 km, is so short in duration of time that favourable weather 
conditions almost simultaneously at all the sites might be just as a stringent require- 
ment as in the case of the geometric mode. (Short arc mode using RF radar 
may alleviate the weather dependency but is negatively compensated by more serious 
refractional problems and more complex active grid stations). 

7. Conclusions and Applications . -- Ranging with 10 cm and 500 observa- 

tions per station can recover relative positions well (Cj... = 4 cm and Ivj... j< 3cm). 
Unit efficiency can be achieved with fewer obse^:^wations (50-100). Expected 

improvements in the^^ ranging accuracy (to 1-2 cm) and in the corresponding pre- 
cision makes the proposed system an excellent candidate for geodetic and geo- 
dynamic applications. , As far as the mode of operation is concerned in case of a 
laser system the following trade-offs need to be considered: 

The likeliness of having- more or less favourable weather conditions at 
4 distant sites in case of the short arc mode (possibly with a single 
satellite and non simultaneous ranging) vs. the feasibility to overcome 
instrumental problems in the geometric mode (airplane and satellite with 
simultaneous ranging). 
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In case of an RF system neither of these problems are critical, and 
the decisive factor is whether systematic errors effecting the RF systems 
can be reduced to the level of those effecting the laser systems. 

Possible candidates as users of a Close Grid Geodynamic Measurement 
System (CLOGEOS) are: Solid Earth - motions near plate boundaries, subsidence 

and uplift, regional strain measurements, horizontal motions, and dilatancy near 
faults, post earthquake resurvey, regional tidal loading, volcanism associated 
motions, surface motions on unstable slopes, geodetic surveys. Cold Regions - 
dynamics of pack ice and ice islands snow /ice motions in major ice sheets, profile 
and flow of glaciers, surface motions in permafrost. Marine Geodesy - positioning 
of ocean bottom geodetic reference frame, positioning or tracking of surface bouys. 
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Attachment 3 


USSR Academy of Sciences 
Astronomical Council 

International Seminar "New Methods of Space Geodesy" 
Leningrad, 24 - 30 November, 1975 


PRELIMINARY PROGRAM 
25 November 


Morning Session (10 am - 2 pm) 

The Opening Ceremony - Balcratov - Greetings 

A.H. MASE’iTCH, N.P. ERPYLYOV, S.K. TATEVYAN (USSR) - Review of Space 
Geodesy Programs as Realized by the Astronomical Council of the USSR 
Academy of Sciences in 1970-1975. 

J. KOVALEVSKY (France) - Installation du Centre d' Etudes et de Recherches 
Geodynamiques et Astronomiques (CERGA) . 

T.J. KUKKAMAKI (Finland) - Utilization of the 890km Long Geodimeter 
Ti'averse in Space Geodesy. 


E.P. FYODOROV (USSR) - On the Observational Methods used for the Earth 
Rotation Studies . 

EITSCHBERGER (FRG) - On Problems of Accuracy of World Geodetic Data from 
Satellites . 

W. PACHELSKI (Poland) - Results of the Analysis of Laser Ranging Measurements 
.and Synchronous Photographic Observations of GEOS-B (1968) b 5 r the 
Successive Adjustments Method. 

Afternoon Session (4 pm - 6 pm.) 

Y. KOZAI (Japan) - Orbital Elements of GEOS-A and -B by Use of Laser 
Observations . 

A. DINESCU, N.RADULESCU (Rumania) - Sur -la Detei'mination preliminaire 
de la station de Bucarest dans le systeme "The Standard Earth. " 
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KISSILEV and BIKOV (USSR) - Orbital Elements from Direct Observations . 

M, BURSA (CSSR) - The Satellite Altimetry and the Scale Factor of the 
Geopotential. 

F.NOUEL (France) -Traitement des mesures et resultats de Geodesic Spatiale 
par recepteur Doppler . 


26 November 
Morning Session (10 am - 2 pm) 

YU.L. KOKURIN, V.K. ABALAKIN (USSR) - On Potentialities and Some 
Results of the Laser Ranging to the Moon, 

I. 1. MUELLER (USA) - Aspects of Positioning using Satellite Borne Lasers. 

V.V. ZLOTIN (USSR) - On Necessary and Practicable Accuracy of Accounting 
for the Light Velocity Variation in the Atmosphere in Laser Ranging to 
AES and to the Moon. 

G. KARSKY (CSSR) - On the Problem of Reduction of Heterogeneous Satellite 

Observations in a Local Network. 

J. KOSTELECKY (CSSR) - Problems of Accurate Reduction of Observations to 

Synchronous Time Moments . 

J. KAKKURI (Finland) - The Finnish Stellar Triangulation Net as a Geodetic 
Control for tlie First Order Terrestrial Trianguiation. 

J. KABELAC (CSSR) - First Realizations of the Trianguiation Project 
using High-altitude Targets in the CSSR. 

M. V. PAUNONEN, A.B. SHARMA (Finland) - Satellite Laser Transmitter and 

Receiver, Technical Solution and Test Results . 

Afternoon Session (4pm - 6 pm) 

H. KAUTZLEBEN, CL. ELSTNER, G. HEMhlLEB, H. MONTAG (GDR) - 

Complex Studies in the Planetary Dynamics of the Earth. 

N, CAPITAINE, L. SAINT GRIT (France) - Variations de la latitude et 

longitude de la station Doppler du CERGA, 

■ N.L. MAKARENKO (USSR) - On the Accuracy of Geometrical Satellite Method 
as Applied to Constructing a Regional Geodetic Network, 
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Morning Session (10 am - 2 pm) 

L.P. PELLINEN, O.M. OSTACH, G.V. DEMYANOV (USSR) ~ On Prospects 
of using the combined Satellite, Gravimetric and Astrogeodetic Data for 
Determination of the Figure and the Gravity Field of the Earth and their 
Time Variations 

L. R. KOGAN, V.I. KOSTENTCO, L.I. MATVEYENKO (USSR) - On Potentialities 

of the Radio-intei-ferometric Facility of the Institute for Space Research as 
Applied to Geodesy and Astronomy. 

W.H. CANNON, R.B. LANGLEY, W.T. PETRACHENKO, N.W. BROTEN, 

D.L. FORT, T.H. LEGGE (Canada), P.A. BARBER, M.J. QUIGLEY 
(England) - Geodetic and Astronometric Measurements using the 
Algonquin- Chilbolton Long Baseline Interferometer. 

P.E. ELIASBERG (USSR) - On Interfering Parameters Affecting the Solution 
of the Problem of Combined Determination of the Earth Figure and Gravity 
Field . 

V.S. GUBANOV, YU. S. STEELETSICY, N.D. UMARBAYEV, B.A. FIGARO 
(USSR) - On Prospects of Solution of Astrometric Fundamental Problems 
by use of the VLSI and Special Space Experiments . 

M. L. LIDOV, YU. F. GORDEYEVA (USSR) - On the Mascons’ Influence on 

Determination of the Moon's Gravitation Coefficients. 

HALMOS, F., ADAM, J., ALMAR, L, FEJES, I. (Hungary) - An Application 
of Radioteclmic Methods (Doppler Measurements) to AES Observations 
for Solution of the Geometrical and Dynamical Problems of Space Geodesy. 

G. BALMINO, B. MOYNOT (France), CH. REIGBER (FRG) - Models de 
potentiel terrestre GRIMl. 


28 November 
Morning Session (10am - 2 pm) 

M. S. PETROVSKAYA (USSR) - On Construction of the Everywhere Convergent 

Geopotential Expansion. 

V. V. BROVAR (USSR) - Coordination of the Satellite and Gravimetric Obser- 
vations in Calculations of Harmonic Coefficients of the Potential of the 
Ellipsoidal Earth. 

GH. VASS (Rumania) - On the APU-75 Network Adjustment. 

N. GEORGIEV, B . SHUSTOV (USSR) - On Mechanical Approximation of AES 

Orbits by Means of Power Series. 
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M. SOIAEIC (Yugoslavia) - On Determining the Distance between Two 
Terrestrial Surface Points by Use of Two AES . 

A. CAZENAVE (France) - Determination des Coefficients des inarees 
oceaniques a partir d'observation des satellites. 

L.K. LAUCENIEKS (USSR) - The General Theory of the One-parameter 
Mobile Barrier. 

IVANOVA (USSR) - Improvement of Orbit using Photographic Obseiwations of 
Planets . 


Closing of Seminar 

Proceedings of the seminars will be published by the Astronomical 
Council, as Vol. XV of the Observations of Artificial Satellites. 
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